Anyonic braiding phases in a rotating strongly correlated photon gas 
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We present a theoretical study of a rotating trapped photon gas where a Laguerre-Gauss laser 
pump with a non-zero orbital angular momentum is used to inject rotating photons into a cavity with 
strong optical nonlinearity. The Laughlin-like few-photon eigenstates appear as sharp resonances in 
the transmission spectra. Using additional localized repulsive potentials, quasi-holes can be created 
in the quantum Hall liquid of photons and then braided around in space: an unambiguous signature 
of the many-body Berry phase under exchange of two quasi-holes is observed as a spectral shift of 
the corresponding transmission resonance. 



Quasi-particles with fractional statistics are among the 
most fascinating discoveries of contemporary condensed- 
matter physics [l|, and are raising an ever growing 
excitement in view of topological quantum computation 
applications 0. In the last decade, nonlinear optical 
systems are emerging as an outstanding newplatform 
to study many-body physics in fluids of light [J] : super- 
fluid hydrodynamic effects have been investigated with 
unprecedented detail @ and strongly correlated states 
like Mott insulator @ or Tonks- Girardeau gases 0, 13 
are under active experimental study. In the present work, 
we present a theoretical study of a rotating trapped pho- 
ton gas: A Laguerre-Gauss laser pump with a non-zero 
orbital angular momentum is used to inject rotating pho- 
tons into a cavity with strong optical nonlinearity, and 
the Laughlin-like few-photon eigenstates 0, [l^l appear 
as sharp resonances in the transmission spectra. Addi- 
tional localized repulsive potentials are used to create 
quasi-holes in the quantum Hall liquid of photons and 
then to braid them around: an unambiguous signature 
of the many-body Berry phase pH under exchange of 
two quasi-holes is observed as a spectral shift of the cor- 
responding transmission resonance. 

The close analogy between the Coriolis force in a ro- 
tating reference frame and the Lorentz force on a charged 
particle in a magnetic field is well-known from textbook 
mechanics and in recent years has been widely exploited 
to study the analogs of the Meissner effect in superflu- 
ids of neutral atoms and the appearance of quantized 
vortices at higher rotation speeds Extension of 

this research into the quantum Hall regime is actively in 
progress using synthetic gauge fields stemming from the 
Berry phase accumulated by an optically dressed atom 
while adiabatically moving in space [Tsj . 

In the photonic context, even though optical vortices 
in nonlinear optical media have received a great atten- 
tion since the earliest works in fluids of light , the in- 
terplay of the orbital angular momentum of light with 
strong photon-photon interactions at the single quan- 
tum level has not been investigated yet. Here, we ex- 
plore this physics in a single cavity geometry: resonant 



transmission peaks are associated [1, [T^ to strongly cor- 
related few-photon states having an excellent overlap 
with Laughlin states [isj . Even more remarkably, the 
many-body Berry phase connected to the braiding of 
one quasi-hole around another [l6j turns out to be di- 
rectly observable as a spectral shift of a resonant trans- 
mission peak, hopefully without the interpretation diffi- 
culties arising from competing effects in electron gas ex- 
periments based on interferometry P, [T7| . Furthermore, 
in contrast to previous studies of quantum Hall states 
of photons [13, [H, nil, ou'^ proposal does not require 
sophisticated photonic configurations and the synthetic 
magnetic field is generated by the non-zero angular mo- 
mentum of the pump. 

The physical system we are considering is a single op- 
tical cavity with cylindrical symmetry consisting of a 
pair of spherical mirrors and containing a slab of non- 
linear medium as sketched in Fig. [IJa). Transverse 
modes with a given longitudinal mode number Afz along 
the cavity axis z can be described as the eigenstates of 
an isotropic two-dimensional harmonic oscillator of fre- 
quency uj ~ y^2c^/RL, L being the central distance be- 
tween the two mirrors and R their radius of curvature. 
Neglecting for simplicity light polarization effects and 
restricting ourselves to a single longitudinal mode, the 
many-body dynamics of cavity photons can be described 
in second quantization terms via the field Hamiltonian Q 
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ftF(r,t)*^(r) +ni^*(r,t)^(r)| . (1) 

Photons in the given longitudinal mode are described by 
the two-dimensional quantum field 'l'(r) satisfying two- 
dimensional bosonic commutation rules [^^(r), >I'^^(r')] = 
^(2)('j. _ 'pj^g confinement between the two mirrors 
is responsible for the finite photon rest frequency ujc = 
ckMzI L and its mass ruph = hujcl <? and the mirror cur- 
vature provides the harmonic trapping potential of fre- 
quency w. Of course, this same Hamiltonian can be used 
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to describe a variety of other configurations, e.g. solid- 
state planar microcavities with a suitable lateral pattern- 
ing [20, , or even hybrid set-ups with a spherical fiber- 
tip mirror facing a planar DBR mirror |22| . The addi- 
tional potential Vqh (r, t) will be taken as a sum of repul- 
sive delta-shaped potentials centered at time-dependent 
positions rf)(t), F,,(r,t) = E.=i K '5(2)(r - rf ^(t)) 
and will serve to create the Nqh quasi-holes in the fluid: 
thanks to the fast modulation speed of optical beams, all- 
optical techniques as the ones demonstrated in [l^, [23] 
can be used to generate these time-dependent potentials. 
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FIG. 1: (a) Sketch of the experimental setup, (b) Steady-state 
relative probability of two-particle excitation P2/A (solid 
line) and of three-particle excitation P3/P2 (dashed line) in 
the absence of repulsive potentials under a monochromatic 
Laguerre-Gauss pump LGq and LGq, respectively, (c) Over- 
lap of the two- (solid line) and three-photon (dashed line) 
amplitudes with the corresponding Laughlin wave function. 
System parameters: Sl/oj = 1, Qni/Puj — 4, -y/to — 0.01, 
IF/'y — 0.1. The three-body calculation is performed within 
the LLL approximation. 



The strength of photon-photon interactions is quanti- 
fied by the interaction coefficient gni proportional to the 
X^^"^ nonlincarity of the medium embedded in the cavity: 
in order to maximize its value, most promising choices 
may be a cloud of optically dressed atoms in a Ryd- 
berg EIT configuration [2^ or, in a solid-state context, 
a quantum well layer with an excitonic optical transition 
strongly coupled to the cavity photon [^[2^ [26j. 

Injection of photons into the cavity by a coherent pump 
laser is described by the last line of the Hamiltonian, the 
spatio-temporal profile of the laser being fixed by the 
function F{r,t): in the following, we shall restrict our 
attention to the case of monochromatic pumps at fre- 
quency ujp^^ with the spatial profile of a Laguerre-Gauss 
beam LG™ with orbital angular momentum m centered 
on the cavity axis. Photon losses at a rate 7 are included 
as Lindblad terms in the master equation for the density 



matrix [27|. 

In the following, it will be useful to describe the sys- 
tem from a reference frame rotating at angular frequency 

around z. To this purpose, one can either add to 
the Hamiltonian ([T]) a single term proportional to the 
total angular momentum Lz along the rotation axis z, 
Hn — H—^Lz, or include the centrifugal force as a reduc- 
tion of the effective trapping frequency ^ — il^ and 
then separately account for the Coriolis force in terms of 
a vector potential A(r) = mp^ fiz x r minimally cou- 
pled to the photon momentum as —ihV — ?> —ihV — A(r). 
While the former formulation is most convenient in cal- 
culations, this latter one emphasizes the close analogy 
with the dynamics of a charged particle in a magnetic 
field. Of course, in the rotating reference frame, the po- 
sitions Ti of the delta potentials have to be accordingly 
rotated back by an angle ili with respect to the labora- 
tory frame ones and the pump frequency is shifted 
to LOp = ujp^^ — mfl. In the following, pump frequen- 
cies will be measured from the photon rest frequency as 
Aujp = ujp — ujc- 

Based on the Hamiltonian ([T]), we now discuss how 
it is possible to generate a Laughlin state of photons in 
the cavity without repulsive delta potentials (Vo = 0). 
We start by considering the isolated system Hamiltonian 
"Ho in the frame rotating at D, in the absence of driv- 
ing (F = 0) and losses (7 = 0). When ^ uj, this 
Hamiltonian is seen to be formally identical to the one 
describing the fractional quantum Hall (FQH) physics of 
interacting electrons in a magnetic field, if one replaces 
the Coulomb interactions with the present contact inter- 
actions. In particular, the A'^-particle ground state of Hci 
is well represented by the = 1/2 bosonic Laughlin wave 
function 



^FQHizi,...,ZN) 



(2) 



where Zj — (xj + iyj ) / £ arc the complex particle coor- 



dinates in units of the oscillator length i — ^Jh/mphto 
[lH . This wave function is composed of lowest Landau 
level (LLL) wave functions, has a total angular momen- 
tum M = hN{N — 1) and is separated from the excited 
states by an energy gap approximately given by / Airt"^ 
in the low gni limit where the LLL approximation is valid 
(cfr. Sec. I of the Supplemental Material). 

Inspired from our previous work [l3|, we propose to 
create such a Laughlin state optically by exploiting the 
driven-dissipative nature of the photonic system. A 
Monte Carlo wave function calculation of the steady-state 
density matrix [13, shows that the cavity transmis- 
sion under continuous- wave pumping and losses displays 
a peculiar dependence on the pump frequency Awp: for 
sufficiently low photon losses, sharp resonance peaks cor- 
responding to the eigenstates of the isolated system ap- 
pear. An effective choice to excite the Laughlin state is 
to use a Laguerre-Gauss LG™ incident beam with orbital 
angular momentum m = A" — 1 chosen to match the value 
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A'' — 1 of the angular momentum per particle of the target 
A/'-particle Laughlin state. The efBciency of this strategy 
is illustrated in Fig. [ljb,c) for the N ~ 2,3 cases. 

In panel (b), we show a simulated spectrum of the 
relative probability Pn/Pn-i of having N photons in 
the cavity as a function of pump frequency Awp in the 
frame rotating at = w. Plotting the ratio Pn/Pn-i 
allows to put in better evidence the role of the A^- 
photon states over the background of lower states: As 
usual , the position of the transmission peak is related 
to the A'^-body cigenenergies by the resonance condition 
= ijji^)/N = ujc+ijJ. For both A^ = 2,3, the main 
peaks at /S.ujp/uj = 1 correspond to an A^-photon tran- 
sition from vacuum to the lowest A^-particle eigenstate 
of at energy hu)^^^ = Nh{ujc + uj) fixed by the zero- 
point motion in the harmonic potential. This suggests 
that particles are non-overlapping in this state. 

As a check of the Laughlin nature of the 
state, we can look at the overlap C'(^'_fq//, $) = 
|(*fqh|^')|^/(*fqh|*fqh)($|'5) between the A^-photon 
amplitude $(zi, . . . , zjv) ~ Tr[pss>I'(zi) . . . ^(zat)] and 
the target Laughlin wave function ^fqh- As we dis- 
cussed in in the weak driving limit £F/j <^ 1, the 
A^-photon amplitude gives in fact the many-body wave 
function of the single A^-particle pure state reached 
by the system and is experimentally accessible from 
repeated measurements of the field quadratures of the 
transmitted light. Its dependence on the pump frequency 
Aujp is shown in Fig. [Ijc): as expected, the maximum 
overlap is obtained at Aujp/uj = 1; its peak value larger 
than 99.5% confirms that the generated state is basically 
the A^-particle Laughlin state. As angular momentum 
of the rotating gas is continuously replenished by the 
pump beam, the photon system is much less sensitive 
to trap anisotropics than atomic clouds: as a result, the 
overlap with the Laughlin state is still w 97% for trap 
anisotropics as large as {ujx ~ ^y)/{^x + ^y) ~ 0.02. 

We can now turn to the generation of quasi-holes in 
the quantum Hall liquid. This can be done by adding 
localized repulsive potentials to pierce holes in the photon 
gas. As indicated in Eq. dT]), their strength is denoted 
by Vo- Their position r-^'(t) is assumed to be rotating 
at an angular frequency in the laboratory frame, in 
order to be stationary at in the frame rotating at f2. 
In the absence of pumping and losses, and for — oj, 
the ground state of the one-quasi-hole Hamiltonian T-Lq 
is successfully represented by the single quasi-hole wave 
function 0,1, [11. 

'^o{zi, ■ ■ .,Zn) = Y\.(^i " ^FQh{zi, ■ . . ,ZAr), (3) 
i 

where Zo = roe^^°/£ is the complex coordinate of the 
quasi-hole (cfr. Sec. I of the Supplemental Material). 
Another quasi-hole sitting, e.g., at the center of the trap 
Zoo = can be included via a second delta-function po- 
tential in the two quasi-hole Hamiltonian Ti^f . Again, in 
the fl ^ UJ limit, the ground state wave function can be 



written in the simple form 

■^oo{zi,...,Zn) = - Zo)zi'ilFQH{zi,...,ZN)- (4) 

i 

The crucial point of our proposal is to relate the braid- 
ing phase observed in the reference frame rotating at the 
trap frequency uj to the time-independent energy spec- 
trum in the frame rotating at slightly lower n = uj — 6fl. 
In the frame rotating at w, the quasi-hole at Zo is in 
fact slowly rotating at frequency Sfl in the backwards 
direction: provided SD, is small enough, this process is 
equivalent to adiabatically looping the quasi-hole at Zo 
along the circle of radius Tq following the position of the 
localized potential. As a result, after a rotation period 
T = 2Tr/SD,, the quasi-hole will return to its original po- 
sition |Tl| , with the single (double) quasi-hole wave func- 
tion (*I'oo) having acquired a Berry phase 0^ {4'b') in 
addition to the trivial dynamical phase E^T {E^°T). 

When observed from the reference frame rotating at 
where the localized potentials are fixed in space, the time 
evolution reduces for any t to the phase t. At time 
t ~ T when a rotation is complete, the wave functions in 
the two frames have to coincide again, which establishes a 
relation between the energy difference AE°'°° = E°'°° — 
E^°° and the many-body Berry phases (f>°g°° , 

A 77^0.00 

0r°=27r-^ (mod27r). (5) 

This relation holds for both quasi-hole states in a quan- 
tum Hall liquid, as well as in a non-interacting system 
(cfr. Sec. HI of the Supplemental Material). As it relates 
the many-body Berry phase to spectroscopically observ- 
able quantities such as the energies, it will be the basis of 
the measurement scheme we are now going to illustrate. 

As a first step, we wish to numerically confirm Eq. 
([5]) for the isolated system. To this purpose, we look for 
the ground state wave functions in the rotating frame 
at n where the quasi-holes are fixed in space and the 
Hamiltonians and H^f are time- independent. Their 
overlap with the analytic wave-functions ([3]) and (|4]) for 
ri/w = 0.99 is shown in Fig. [2][a) as a function of the 
position To of the exterior quasi-hole. In the lower panel, 
we show the value of the braiding phase ^bi- = 4>°b ~ '■ 
as introduced in [l^, this is the difference between the 
many-body Berry phases acquired by the ground state 
wave functions in the presence of single and double delta- 
function potentials. In this panel, the value of 4>Br ex- 
tracted via Eq. ([S]) from experimentally accessible quan- 
tities is compared with the result of a direct calculation 
of the Berry phases from the analytical wave functions 
(O and (gl) (cfr. Sec. II of the Supplemental Material). 
The agreement is excellent up to a radius ~ i, i.e. 
when the quasi-hole potential starts exiting the cloud: at 
this point, the repulsive delta potential is no longer able 
to sustain the quasi-hole state and the overlap shown in 
panel (a) suddenly drops. 

The optical generation of the quasi-hole states is il- 
lustrated in Fig. 131 a monochromatic Laguerrc-Gauss 
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FIG. 2: (a) Stars * (circles o) show the overlap of the N = 2 
lowest energy eigenstate of Hq {'H'q) with the one (two) 
quasi-hole wave function ^l/o (^Poo) as a function of the po- 
sition ro of the exterior quasi-hole. (b) Estimation (•) of the 
braiding phase calculated via Eq. ([5]) compared with the an- 
alytical result from the quasi-hole wave functions ^l/o, ^'oo 
(solid line). System parameters g„i/l^uj = 4, Sl/oj — 0.99 and 

Vo/f hu = 100. 
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FIG. 3: Steady-state relative probability P2/P1 of two- 
particle excitation for "Hq with a Laguerre-Gauss LGq pump 
(solid line), and for Hn" and an LGq pump (dashed line). Ver- 
tical solid (dashed) lines correspond to half the eigenenergies 
of the isolated system. System and pump parameters as in 
Fig. [H with Q./u] = 0.99, To = 0.4 and Vo/l^fiw = 100. 



beam is shined on the cavity in the presence of the re- 
pulsive potentials rotating at a frequency Vt/uj = 0.99. 
An efRcient choice for the orbital angular momentum of 
the pump is to use the closest integer to the angular 
momentum per particle of the target state. The steady- 
state density matrix is numerically calculated via Monte 
Carlo wave function technique in the frame rotating at 
51 where the Hamiltonian is time-independent. The dif- 
ferent curves in the figure show the spectrum of P2/P1 
as a function of pump frequency Aw^: the solid (dashed) 
curve refers to the one (two) quasi-hole Hamiltonian 
{TLq) including one (two) localized repulsive potential. 
The most relevant features are the lowest frequency peaks 
at Auip/uj = 1.0195 on the solid line and 1.0297 on the 
dashed one: their identification with quasi-hole states is 
confirmed by the excellent overlap w 99% of the two- 
photon amplitude with the analytical wave functions in 
([3]) and (H]). Remembering that these two-photon peaks 
are located at half the energy E'^2°° oi the two-photon 
eigenstate, it is then straightforward to extract via ([5]) 
the value of the many-body Berry phase when one quasi- 
hole at To is braided around another one located at the 
center of the trap: as one can see from the small circle in 
Fig. IH^b), the accuracy of this simulated measurement 
is excellent. An analogous calculation of the position 
of three-photon peaks (performed within the LLL ap- 
proximation) suggests the efficiency of the measurement 
protocol also for states with larger number of particles 
where the theory of the fractional quantum Hall effect 
predicts a constant braiding phase of tt (cfr. Sec. IV of 
the Supplementary Material) 0, [l^ . 

In conclusion, we have proposed and characterized an 
all-optical scheme to generate and manipulate quantum 
Hall liquids of strongly interacting photons in a nonlinear 
optical cavity. Quasi-holes in the fluid can be pierced and 
braided with repulsive potentials and the corresponding 
many-body Berry phase can be detected from the spec- 
tral shifts of the resonant transmission peaks. Extension 
of this work to more complex configurations involving e.g. 
light polarization degrees of freedom may open the way to 
observe anyonic excitations with non-Abclian statistics. 

We arc grateful to A. Imamoglu, R. Santachiara and T. 
Volz for stimulating discussions. The authors acknowl- 
edge financial support from ERC through the QGBE 
grant. 
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I. DEPENDENCE OF EIGENSTATES AND 
EIGENENERGIES ON THE RELATIVE 
ANGULAR FREQUENCY SQ 

In this section we justify the adiabatic following of the 
ground state when the quasi-hole potential is braided 
around in space. 

Let us start from the Hamiltonian with no re- 
pulsive localized potentials. The energies of the eigen- 
states of Hii are plotted in Fig. [TJa) as a function of 
(50. For Jfi — i> 0, we have a sequence of massively de- 
generate states: the lowest manifold corresponds to wave 
functions obtained as the product of a Laughlin state 
^ FQHi as defined in Eq. (3) of the main text, times any 
polynomial symmetric under zi ^ Z2 and is separated 
from the higher manifold by a gap mostly determined by 
the interaction strength gni- Within the lowest Landau 
Level (LLL) approximation valid for small gni/^"^ ^ 
the gap is proportional to gni and close to the integral 
9ni j \iPq{v)\'^(Py/2 = Qnil-iTrP, V'o being the harmonic 
oscillator ground state wave function. For finite 5^1 , the 
states of the first manifold are split according to their 
total angular momentum. The non-degenerate ground 
state has the lowest angular momentum hN{N — 1) and 
is given by the Laughlin state '^fqh- The first excited 
state has one more unit of angular momentum and is 
separated from the ground state by h5^l. 

We have then performed a numerical diagonalization of 
the Hamiltonian Ti!l{°° in the presence of one and two lo- 
calized repulsive potentials piercing quasi-holes through 
the quantum Hall liquid. The result for the two quasi- 
hole case is shown in Fig. [ijb): for small but non-zero 5Vl, 
the ground state is separated from the first excited state 
by a small amount close to hSVl. The next manifold lies 
higher in energy by an amount determined by gni- States 
involving non-lowest-Landau-level components are found 
at even higher energies of the order of fiw. The quasi- 
hole nature of the ground state is numerically confirmed 
by a calculation of its overlap with the two quasi-hole 
wave function, as given by Eq. (4) of the main text. The 
higher states of the manifold correspond to wave func- 
tions obtained by multiplying this two quasi-hole wave 
function by arbitrary symmetric polynomials. 

The fact that the two quasi-hole state is the ground 
state even for non-zero values of 5Vl confirms that holes 
are indeed braided around the quantum Hall fluid with- 
out creating excitations provided braiding is performed 
slowly enough. As we have seen in the main text, the 
overlap with the quasi-hole wave function suddenly drops 
at some value of the quasi-hole position Tq . An analogous 
phenomenon is visible here as a function of angular fre- 
quency 5VL in Fig. [IJc): overlap suddenly drops at the 
first avoided crossing point between the ground state and 
a lower angular momentum state from the higher mani- 
fold. 
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FIG. 1: (a) Energy of the lowest N = 2 particle eigenstates 
of Hn as a function of relative rotation frequency S^l = lu — 
Q in the absence of the repulsive quasi-hole potentials. For 
clarity, only the states of total angular momentum Lz < 5 are 
displayed, (b) Energy of the lowest N = 2 particle eigenstates 
of the two quasi-hole Hamiltonian H'q as a function of relative 
rotation frequency Si} — uj — il. For clarity, only states at 
low average angular momentum (Lz) < 8 are displayed, (c) 
Overlap of the ground state of Hq' with the analytical two 
quasi-hole wave function, Eq. (4) of the main text. System 
parameters: g„i/£'^uj — 4, Vo/I!'^ fiu! — 100. The exterior quasi- 
hole is positioned at ro = 0.4^". 



II. BERRY PHASE FOR ONE AND TWO 
QUASI-HOLE STATES 

Here we prove a simple relation between the expected 
value of total angular momentum (L^) and the Berry 
phase (jjB acquired by a many-body wave function lying 
in the LLL, when a quasi-hole is looped along a circle 
with or without a second quasi-hole sitting at the center 
of the circle. 

The generic form of the A^-particle wave function we 
consider is 



^ ^ '^ni,...,nN ^1 ' 
ni ,. . . ,njv 



^{zi,...,zn;zo) 

X (zi-0o)...(zAr-Zo)e-2:;=il-.lV2^ (1) 
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where Zj = {xj + iyj)/^ are the complex particle coor- 
dinates, Zo is the quasi-hole coordinate and tij are non- 
negative integers. The only assumption we are making is 
that X^jLi 'T-i = ^ is constant, meaning that the part of 
the wave function shown in the first line is a homogenous 
multivariate polynomial of degree A. A second quasi-hole 
term zi . . . z^v can be included into this part by letting 
each Tij increase by 1. 

For simplicity, from now on we will omit displaying 
the exponential factor which is common to all terms in 
the summation over Uj and use the shorthand notations 
{z} = zi,...,ZAr, and {n\ = ni,. 
expand the product (zi — Zo) . . . (zjv 
Zo) ■ ■ . (zN — Zo) — J2i^QdiZo. Upon these changes, the 
wave function becomes: 



,nN- We will also 
Zo) as follows (zi — 



While taking the derivatives dg^ with respect to the angu- 
lar variables (we have set Zj = \zj \ e*^^ ), we can use the 
chain rule to isolate the contribution of the "^^^gdiz^ 
term. In this way, the second line of Eq. Q becomes 
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Let us now calculate (f>B for a quasi-hole at Zo — roC 
with fo = To/l fixed: 



^B = -^J^ J d{z}d{z*}-^*dg^^d0o = 

= 2njd{z}d{z*} ^ c^„,jC{„|zr^zr^ 



/£ To get this result one simply has to recall that d 
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where c?{z} = dzi . . .dz^ and dg^ = d/dOo is the partial 
derivative with respect to 9o- Now let us look at the 
expected value of total angular momentum: 
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with (s 



+ sn) = N - I. Inserting ^ 



into Eq. and reordering terms we obtain 



{"},{»'} 
i.i' 
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X {-lh){dg,+ . . . +de,)Y, C{n}zT- ■ • ^'^o(4) 
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In order for the integral not to vanish, the total power 
of conjugated coordinates should match that of unconju- 
gated ones: A + N-l' = A + N-l, implying I' ^ I. As 
a result, one gets 



- 27ri^ = -2n{A + N)Y, c^„,jC{„j d{z}d{z*}zf\ . . z"^ 

+ cl^,yCi^Jd{z}d{z*}zr\..z-^Y^\d,\^^o- (7) 

I 

The first term apart from the factor —2tt{A + N) is the Thus, finally we have 
normalization of the wave function ([2]), which wc take to > 

be 1. The second term is 0b we have found in Eq. (ps = -211^^ + 2t:{A + N). (8) 

h 

This result can be explicitly checked for the simplest 



3 



N = 1 and N = 2 cases, for which one can analytically 
calculate (pB and (L^) without much difficulty. Here we 
do not give the details of the elementary verification and 
only quote some final results for one and two quasi-hole 
states. 

For the case of = 1, the Berry phase for the single- 
quasi-hole wave function oc (zi — Zo)e~'^^' (here ^ = 0) 
is 



a1-° 
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(9) 



while for the double-quasi- hole wave function oc z\ {z\ 

Zo)e"l'=il'/^ (here A = 1) it is 



(10) 
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For the case of iV = 2, the Berry phase for the j/ = 1/2 
Laughlin wave function with a single quasi-hole oc (zi — 

Z2f{zx - Zo){z2 - Zo)e-(l"il'+l"^l')/2 (here ^ = 2) is 
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while for the double-quasi-hole wave function oc (zi — 

Z2)2ziZ2(zi - Zo)(z2 - Zo )e-(l^i l')/^ (hcrc A = 4) it 
is 
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18fg + 14f^ 
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(12) 



To summarize, the relation given in Eq. ([5]) provides a 
direct way to efficiently calculate the Berry phase as soon 
as the coefficients of each term in the wave function are 
known [1]: in the last section we shall see an example of 
its application to states with a larger number of particles 
N up to 6. 

In addition to this, this relation suggests an alterna- 
tive experimental method to determine the Berry phase 
for setups where a precise measurement of total angular 
momentum might be feasible. 



III. RESULTS FOR THE NON-INTERACTING 
SYSTEM 



FIG. 2: (a) Stars * (circles o) show the overlap of the lowest 
energy state wave function of a single particle in the presence 
of one (two) delta-function potential with the one (two) quasi- 
hole wave function as a function of the exterior quasi-hole 
coordinate ro. System parameters: Q,/uj = 0.99, Vo/i'^huj = 
100. (b) Braiding phase calculated from Eq. (|13p is shown by 
dots (•); the curve is the analytical result for (/!>]^° — 0^°° as 
given by Eqs. Q and ((101 . 



The many-body wave function for a larger number N of 
non-interacting particles is in fact equal to the product 
of = 1 wave functions, so that both the energy and 
the braiding phase are just multiplied by a factor N. 

The braiding phase as a function of the exterior quasi- 
hole coordinate Tq is plotted in Fig. [2{a) and compared 
with the difference 0]^° — (j>g°° calculated from the ex- 
plicit expressions in Eqs. ^ and ((TO)) : once again, the 
agreement is excellent as long as the repulsive potential 
is well inside the cloud, Tq < 0.5^. This explanation is 
validated in panel (b) where we show the overlap of the 
ground state with one and two repulsive potentials with 
the corresponding analytical forms of the one and two 
quasi-hole wave functions. 



In order to get a deeper understanding of the braid- 
ing phase in a quantum Hall liquid, it is interesting to 
study the structure of the ground state and calculate the 
Berry phase also in the different case of a non-interacting 
system with gni — 0. 

With no loss of generality, we can start from the = 1 
single particle case and calculate the braiding phase <pBr 
in the presence of one and two delta-function potentials 
using the relation (5) of the main text 

A po,00 

0^ = 2^-^ (mod 2^). (13) 



The optical response of the gni — system is a purely 
linear one: transmission spectra for a Laguerre-Gauss 
LGq mode with one unit of angular momentum are shown 
in Fig. 13] for the single quasi-hole case. As expected, the 
position of the one photon peak matches the ground state 
energy of the isolated system. The two-photon peak in 
the P2/P1 spectrum occurs at the same frequency, as ex- 
pected from the simple fact that the two-photon ground- 
state energy is twice the single-particle ground-state en- 
ergy. The Berry phase calculated from the shift Auo of 
the peak via Eq. is ~ 0.2807r, in good agreement 

with the analytical prediction w 0.2767r from Eq. ©. 
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FIG. 3: Steady-state probability of single-particle excitation 
Pi (solid line) and the relative probability of two-particle 
excitation P2/P1 (dashed line) in the presence of a single 
delta potential localized at ro = 0.4^ under a monochromatic 
Laguerre-Gauss LGj pump. System and pump parameters: 
Q,/uj = 0.99, -i/uj = 0.01, IFh = 0.1. 



IV. RESULTS WITH MORE THAN TWO 
INTERACTING PARTICLES 



Some results for larger number of interacting particles 
are shown in Figs. |3]and[S] 

To obtain the braiding phase 0Br analytically for N > 
2 more efficiently, we used the relation (|8]). In Fig. Ufa), 
we compare the estimation for 0Br obtained for = 3, 4 
from the numerical ground state via Eq. (1131) with the 
analytical results from the quasi-hole wave functions. As 
compared to the N = 2 case, it can be seen that the re- 
gion of good agreement extends further in ro , which could 
be expected as the cloud size increases with increasing N. 

Panel (b) displays the analytical predictions for (j>Br as 
a function of Tq for up to 6; unfortunately, available 
computing resources prevent us from directly extending 
the calculation to higher N values for which more sophis- 
ticated Monte Carlo techniques should be used. For in- 
creasing TV, the position of maximum (f)Br shifts towards 
larger ro and its value gets close to tt. This trend is quite 
different from that of the non-interacting system exam- 
ined in the previous section, where the braiding phase 
is simply given by N times the single-particle braiding 
phase. 

For much larger number of particles, the characteristic 
incompressibility of the quantum Hall liquid is expected 
to result in a large plateau at (/)Br = tt. In this limit 
0, the Berry phase acquired through looping of a 
quasi-hole can simply be expressed as 2Tr{N), {N) being 
the average number of particles enclosed by the loop. The 
effect of an additional quasi-hole inside the loop would be 
to displace a fraction v of a particle from an otherwise 
spatially constant density, leading to an extra braiding 
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FIG. 4: (a) Braiding phase (pBr calculated for = 3 (dashed 
line) and = 4 (solid line) from the analytical forms of the 
one and two quasi-hole wave functions, Eqs. (3) and (4) of the 
main text. Corresponding (pBr extracted via Eq. (|13p from 
the numerical ground state of Hq"" within the LLL approx- 
imation for A'^ = 3 (•) and N = 4 {<>). System parameters: 
Q./LJ = 0.99, g^i/fuj = 4, and Vo/ftiw = 100. (b) Braid- 
ing phase (j>Br calculated analytically for A'^ = 2, 3, 4, 5, and 6 
from the analytical forms of the one and two quasi-hole wave 
functions, Eqs. (3) and (4) of the main text. Line thickness 
increases with A'. 




FIG. 5: Steady-state relative probability P3/P2 of three- 
particle excitation in the presence of a single delta potential 
localized at ro = 0.4^ under a monochromatic Laguerre-Gauss 
LGo pump. System and pump parameters: = 0.99, 



u = 1, 7/0; = 0.01, IF/'-i = 0.1. 
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phase of I-kv. For the Laughlin and quasi-hole states 
considered in the present work, the fraction is v = 1/2, 
which results in a braiding phase of tt. 

In Fig. [S]we simulate the optical excitation of an = 3 
particle, one quasi-hole state using a Laguerre-Gauss LGg 
pump with three units of angular momentum. At the 
position of the transmission peak at Awp/oj w 1.0297, 
the overlap between the three-photon amplitude and the 
one quasi-hole wave function has the excellent value of 



98.9%. The Berry phase calculated from the frequency 
shift Acjo of the peak via ([T^ is « O.IGStt. This value is 
very close to the analytical prediction « 0.1607r obtained 
from the analytical form of the one quasi-hole wave func- 
tion, Eq. (3) of the main text, in the three-photon case. 
This suggests that our scheme to measure the many-body 
braiding phase can successfully extend to states with a 
larger number of particles. 
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